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Field dependence of magnetic correlations through the polarization flop transition in 
multiferroic TbMnOs : evidence for a magnetic memory efl"ect 
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The field-induced multiferroic transition in TbMnOs has been studied by neutron scattering. 
Apart strong hysteresis, the magnetic transition associated with the flop of electronic polarization 
exhibits a memory effect: after a fleld sweep, TbMnOa does not exhibit the same phase as that 
obtained by zero-fleld cooling. The strong changes in the magnetic excitations across the transition 
perfectly agree with a rotation of the cycloidal spiral plane indicating that the inverse Dzyaloshinski- 
Moriya coupling causes the giant magnetoelectric effect at the fleld-induced transition. The analysis 
of the zone-center magnetic excitations identifles the electromagnon of the multiferroic high-fleld 
phase. 

PACS numbers: 75.30.Ds, 75.47.Lx, 75.40.Gb, 77.80.Fm 



The magnetoelectric effect couples the electric polar- 
ization P to an applied magnetic field H and the magne- 
tization M to an external electric field .E [Tl[2]. Very re- 
cently, giant magnetoelectric coupling has been reported 
for several transition-metal oxides, such as e. g. per- 
ovskite manganites RMn03 with R=Gd, Dy, and Tb 
[H nil spinel chromate CoCr204 [5], spin-chain cuprate 
LiCu202 |6], and huebnerite MnW04 |71 |8]. Com- 
mon to all these systems is a non-collinear magnetic or- 
dering, the key element to understand the multiferroic 
order|9l [10] [11]. Two non-collinear spins Si, Sj at a dis- 
tance Tij induce a spontaneous electric polarization via 
an inverse Dzyaloshinski-Moriya coupling 



P OC Tij X {Si X Sj). 



(1) 



A spiral structure with a finite angle between the spin- 
rotation axis Si X Sj and the propagation vector k , there- 
fore, induces a finite electric polarization, and a large 
magnetoelectric effect may arise from the response of the 
spiral to an applied magnetic field jlOl [12] . 

TbMnOg is particularly well suited for the investiga- 
tion of the oxide multiferroics, as it exhibits a rather 
large electric polarization and as its crystal and mag- 
netic structure is simpler than that of other multifer- 
roics. TbMnOg crystallizes in the orthorhombic symme- 
try, space group Pbnm |T3], and the Mn magnetic mo- 
ments order below Tn=42K in a collinear spin-density 
wave (SDW) with incommensurate (IC) propagation vec- 
tor fcMn ~ (0 0.28 0) |14| [15]. In a second magnetic tran- 
sition at rFE=28K, the Mn magnetic order continuously 
transforms into an elliptical cycloid with the spiral basal 
plane being parallel to the be plane (i. e. the spins ro- 
tate around the third crystallographic axis a) [1^. Ac- 
cording to the general concept of eq. [l] the spiral order 
implies a spontaneous electric polarization P^ along c, 
which is observed experimentally below TpE |3j, render- 
ing TbMnOg multiferroic, see Fig. W\ Applying a mag- 
netic field, giant magnetoelectric effects are observed. A 
critical field Hc=7T along c completely suppresses the 
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FIG. 1: (color online) Sketch of the magnetic spiral structure 
in zero fleld inducing an electric polarization Pc (above). The 
rotation from Pc to Pa in the high-fleld phase is supposed to 
originate in a concomitant flop of the spiral rotation axis from 
a to c (below). 



electric polarization [T71[T8], while a field of 10.5 T (6T) 
parallel to a (b) fiops the electric polarization by 90° from 
Pc to Pa [3j. Since the magnetic modulation in the FE 
high-field phases remains along [010], the 90°-rotation 
of the electric polarization suggests the fiop of the spiral 
rotation axis from a to c, see Fig. IT] [Toj [12] . However, a 
proof that inverse Dzyaloshinski-Moriya coupling accord- 
ing to equation (1) accounts for the giant magnetoelectric 
effects at the field-induced transition is still missing. 

Measuring selected superstructure refiections by neu- 
tron or x-ray diffraction techniques |18l [T9l [20] it has 
been demonstrated, that a magnetic field H parallel to c 
melts the IC ordering and stabilizes a paraelectric phase 
with simple AFM ordering [18]. The fiop of the electric 
polarization for a field along a or 6 coincides with a first 
order transition into a high-field commensurate (HF-C) 
ordered phase with propagation vector q^^ — (0 0.251) 

nil [20]. 

The mechanism in eq. (1) provides for a close cou- 
pling between the dielectric properties and magnetism 
which leads to hybridized phonon-magnon excitations 



[21 [To] [TT] as they were proposed long ago [21]. While 
such electromagnon was proposed on the basis of In- 
frared (IR) optical-spectroscopy on TbMnOa [22], only 
the correspondence of the peaks in the IR spectrum with 
neutron scattering measurements of the IC zone-center 
magnetic excitations [23] documents the mixed phonon- 
magnon character of the excitation. This rather soft elec- 
tromagnon mode possesses the correct symmetry to ro- 
tate the spiral plane from the be to the ah plane and 
appears to be activated when field is applied along a or 
b. 

In this work we follow the antiferromagnetic zone- 
center excitations through the polarization flop transition 
with H applied parallel to a using inelastic neutron scat- 
tering. Our measurements show pronounced changes in 
the magnetic excitations at the critical field, while anal- 
yses of the spectra are consistent with a magnetic spiral 
where the magnetization rotates within the ab plane as 
predicted by the inverse Dzyaloshinski-Moriya coupling. 
Detailed measurements of the propagation vector k as 
a function of magnetic field reveal a magnetic memory 
effect that is driven by domain wall locking. 

Neutron scattering experiments were performed at the 
cold triple-axis spectrometer PANDA (FRM-II, Garch- 
ing) using the same single crystal as in our previous stud- 
ies [23]. The sample was mounted with the scattering 
plane defined by [0 1 0] and [0 01] in a 15 T-cryomagnet 
and the field applied along the vertical a-axis. The 
P6nm-lattice constants are a = 5.302 A, b = 5.857 A, 
and c = 7.402 A, which in the following will be used to 
index all vectors of reciprocal space. Monochromatic neu- 
trons were selected and analyzed using the (00 2)-Bragg 
reflection of pyrolytic Graphite (PG). In the standard 
set-up the energy of the analyzed neutrons was fixed to 
£^/=4.66meV {kf = 1.50 A ), but in order to enhance 
the experimental resolution selected scans were repeated 
with £;/=2.98meV {kf = 1.20A~Y 

Fig. |2] shows the response of the magnetic superstruc- 
ture reflection Q — (Oef,l) to a magnetic field i?a||Q 
for a complete field cycle ffa = OT ^ 12T ^ OT in 
the ferroelectric regime at T = 17 K. The intensity of 
the magnetic reflection increases continuously with in- 
creasing field, see Fig. |2ji. In contrast, up to 8T the 
magnetic modulation described by the zero-field incom- 
mensurability e° = 0.2772(1) remains nearly unchanged. 
With further increase of the field across the critical value 
i/? ~ 9T, £b exhibits a sharp discontinuity, and for 
H > H^ the magnetic intensity is centered at the com- 
mensurate position Q — (0 0.251), see Fig. [2t, and b. In 
the vicinity of iJf, a region of coexistence of low-field in- 
commensurate (LF-IC) and HF-C phases exists, and the 
transition is clearly of first order, in good agreement with 
recent scattering experiments |19l |20] and thermody- 
namic studies [H]. On the downward run a pronounced 
hysteresis is observed, as the system switches back into 
the LF-IC phase at iJ]" ~ 8T [24J. Again, around H^ 
a region of coexistence of both phases is found, which is 
broader than that in the upward stroke. Furthermore, 
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FIG. 2: (color online) (a) Cut through reciprocal space along 
Q = (0 fc 1) for a complete field cycle J/a = 0T-^12T-^0T 
and Ha\\a. Open points mark the positions of the magnetic 
superstructure reflections, dashed lines denote the phase tran- 
sitions between the LF-IC and the HF-C phase as determined 
in Ref. [23] (b) Field dependence of the incommensurability 
Eb for the complete fleld-cycle. The grey shaded regions mark 
regimes with coexistence of both magnetic phases, solid lines 
are included as guides to the eye. (c) Hysteresis of the mag- 
netization measured with iJa||a sensing essentially the align- 
ment of Tb-moments, note that there is a minor hysteresis at 
the LF-IC to HF-C transition. 



upon reentering into the LF-IC phase the magnetic in- 
commensurability £b does not recover its initial value, 
see Fig. |2]3. Switching back from the HF-C phase we 
find £;,(8T|) = 0.2624(3), which is significantly smaller 
than in the upward run, efc(8T|) = 0.2708(1). At zero 



field 



.Oi 



0.2716(1) is smaller than the initial value 
e^ determined before the field sweep. We can only re- 
cover the original zero-field incommensurability by heat- 
ing the system up into the paraelectric SDW-phase and 
subsequent re-cooling it to 17 K. In the FE spiral phase, 
TbMnOg memorizes its magnetic history as the size of 
the incommensurability gjj at zero field depends on the 
magnetic diary of the sample. Measurements of magne- 
tization from a single crystal of TbMnOg with H\\a show 
that there is no large effect that tracks the hysteresis of 
the magnetic wave vector as we found in our neutron 
measurements. The magnetization hysteresis that is as- 
sociate with the transition to the HF-C phase is very 
small, see Fig. |2j;. This indicates that the ferromagnetic 
alignment of Tb-spins that occurs with applied magnetic 
field is of little consequence to the memory effect we find 
here. A more plausible explanation is that the transition 
into the HF-C phase induces domain walls which at least 
partially remain when switching back to the LF-IC phase. 
These domain walls are only suppressed when fully leav- 
ing the spiral phase upon heating [25]. Memory effects 
at combined spin and polarization flop transitions should 
be of general importance in view of future applications 
in magnetoelectric data storage. 

We now turn to the field dependence of the magnetic 
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FIG. 3: (color online) (a) Energy scans at the magnetic zone- 
center Q = {Osb 1) in the LF-IC phase at _ff = OT and in the 
HF-C phase at _H" = 12 T, both recorded with Ef = 4.66 meV. 
The zero-field data are the same as shown in Ref. l23l (b) High- 
resolution data of the 12T-spectrum for Ef = 2.98 meV. (c) 
Raw-data scans for various magnetic fields, and (d) differ- 
ence spectra derived from the data presented in (c). Solid 
lines correspond to fits to the data in (a)-(c), vertical bars 
mark the position of the magnetic excitations at zero field, 
small diamonds label regions with spurious contributions in 
the 12T-spectra. 



excitation spectrum. Inelastic scans at the magnetic 
zone-center Q = (Oefcl) for T = 17K are presented in 
Fig. |3) At zero field, several magnetic signals are de- 
tected. The feature around 4.5 meV corresponds to a 
crystal field (CF) excitation in the Tb-subsystem [23]. In 
the following discussion we will ignore this CF-excitation, 
as up to the maximum field investigated, H = 12 T, no 
significant changes can be observed in the spectra above 
« 4meV, see Fig. |3ji. In the regime of the Mn spin- wave 
excitations below 4 meV the impact of the magnetic field 
is, in contrast, remarkable, see Fig. |3}; and d. 

At zero field, the magnon spectrum consists of three 
different low-energy branches: One mode with very low 
energy, and two modes at finite energies, 1.0 meV and 
2.5 meV, respectively, are observed at the magnetic zone 
center [23]. With increasing field, this characteristic 
structure of the spectrum remains unchanged up to Ha — 
8 T, and the data for 4 T and 8 T can well be described 
using similar parameter sets as those for the zero-field 
data. The observed spectral weight of all modes increases 
following the elastic signal with increasing field, but for 
Ha < 8 T we do not observe significant changes in the 
spin-wave frequencies. Note, that the broad response 
close to uj = OmeV, which has to be ascribed to the low- 
lying spin-wave mode, can not be fully separated from 
the elastic signal. 

With the transition into the HF-C phase the spec- 
trum exhibits prominent changes, see Fig. |3|i. The sig- 
nal around 1.0 meV is completely suppressed, and for 
if" = 12T the spectrum can be decomposed into two in- 
tense excitations centered around 0.5 meV and 2.25 meV, 
and a weaker feature around 3.0 meV. A fourth, rather 
sharp signal around « 1.5 meV, marked by small dia- 



monds in Fig. |3| can be ascribed to a spurious contri- 
bution ^6]. The low-energy spin- wave spectrum of the 
HF-C phase consists of three branches, similar to the 
spectrum in the LF-IC phase |23) . 

The data measured at lOT-data can be described by 
a weighted summation of the 8 T and 12 T spectra, see 
Fig. |3j;, consistent with the co-existence regime of the C 
and IC wave vectors shown in Fig. 2b. 

In the zero-field phase with P\\c longitudinal polariza- 
tion analysis allowed for the unambiguous identification 
of the phason mode and of the Goldstone Boson of the FE 
spiral structure [23] . The applied magnetic field prohibits 
a similar polarization analysis for the P^ HF-C phase. 
However, the systematic survey of the Q-dependence of 
the magnon signal can determine the character of the dif- 
ferent magnon branches in the HF-C phase, as only mag- 
netic polarization perpendicular to the scattering vector, 
S±, contributes to the neutron scattering intensity. In 
Fig. |4| we compare the excitation spectra recorded at 
the three different A-type zone centers Qi = (0 0.251), 
Q2 = (01.751) and Q3 = (0 0.25 3). At these three 
equivalent q-points, the observed magnon frequencies are 
identical, fuvi = 0.44(1) meV, huj2 = 2.20(2) meV and 
huj3 = 3.06(2) meV. However, significant changes can 
be observed in the distribution of spectral weight among 
the three modes. The mode 0^2 is strongest for Q — Qi, 
whereas the branches wi and uj^ are most intense at Q2 
and are remarkably suppressed at Q3. With f{Q) denot- 
ing the magnetic form factor and a denoting the angle 
between Q and the magnetic polarization, the observed 
intensity is given by 



/cx/2(Q)sin2(a). 



(2) 



With the three chosen Q points the geometrical condi- 
tions significantly vary as the angle with c, ac, changes 
: ad — 17.5° at Qi, ac2 — 65.5° at Q2, and Q3 is al- 
most parallel to c, Ucs — 6°. The strong suppression of 
spectral weight for uji upon rotating the scattering vector 
towards c thus implies a magnetic fluctuation polarized 
predominately along c. Indeed, from the observed inten- 

1 : 4.8(4), which is close 
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FIG. 4: (color online) (a) Comparison of raw-data energy 
scans in the HF-C phase for three equivalent Q-positions on 
a logarithmic scale, and (b) on a linear scale. Solid lines 
denote fits to the data and grey vertical bars mark the HF 
magnon-frequencies. Small black diamonds denote regions 
with spurious contributions for Q — (0 0.25 1). 



to the expected value for an entirely c-polarized mode, 
Iq^ : Iq^ = 6.2. The second mode W2 exhibits oppo- 
site behavior. The loss of scattering intensity from Qi to 
Qs is fully explained by the magnetic form factor f{Q), 
implying that uj2 is polarized perpendicular to c. The 
third mode 0J3, although by far the weakest of the three 
magnon excitations, exhibits similar behavior as cui with 
the most intense signal for Q = Q2, see Fig. [ija, and 
we may thus conclude that uj^ also possesses a dominant 
component along c. 

The magnetic structure in zero field is described by a 
magnetic cycloid in the bc-plane [16], and it has been 
shown both theoretically and experimentally that the 
low-energy spin-wave spectrum consists of three branches 
with different character |23l[27| . One of the modes is the 
sliding mode of the spiral and is polarized within the 
spiral plane. The two other modes correspond to mag- 
netic fluctuations perpendicular to the spiral plane and 
are polarized along a. These modes rotate the spiral 
basal plane around c and around b, respectively. The 
latter one is identified as the electromagnon observed in 
IR spectroscopy |22]- This mode thus possesses the sym- 
metry to rotate both the magnetic spiral plane and the 
ferroelectric polarization. 

Applying the arguments of the zero-field analysis to 
the spectrum of the HF-C phase, our observations are 
fully consistent with the field-induced flop of the spiral 
to the ab-plane, as the polarization patterns of the spin- 
wave excitations are rotated by 90°. One of the three 
modes in the HF-C phase, lu2, is polarized within ab, 
and this mode should correspond to the sliding mode of 
the commensurate spiral. Compared to the LF-IC phase 
the phason energy is strongly enhanced, increasing from 
«0.2meV for 7J = OT to 2.20 meV in the HF-phase. The 
large phason energy indicates a strong pinning in the HF- 
C phase resulting in a very anharmonic magnetic modu- 
lation associated with domain boundaries. This further 
corroborates our interpretation that the magnetic mem- 
ory effect arises from domain boundaries introduced in 
the HF-C phase. Furthermore, the magnetic field will ad- 
ditionally trap the phase of the spiral. The two other HF- 
modes are both polarized predominantly along c, i. e. per- 
pendicular to the flopped spiral plane. These branches, 
thus, correspond to the two a-modes of the zero-field spi- 
ral, and following the above argumentation, at least one 
of these modes is expected to couple strongly to an alter- 
nating electric field along c and should be visible in op- 
tical spectroscopy. This mode is the multiferroic electro- 
magnon, but our neutron scattering experiment cannot 
determine which of the two c-modes at high field corre- 
sponds to the electromagnon. We are, unfortunately, not 
aware of an optical study to observe the electromagnetic 
response of TbMnOs in the HF-C phase. 

In conclusion, we have studied the field dependence 
of elastic and inelastic neutron scattering in multiferroic 
TbMn03. The magnetic superstructure reflections ex- 
hibit strong hysteresis across the magnetic transition ac- 
companying the flop of the electronic polarization and, 
most interestingly, there is a clear magnetic memory ef- 



fect. After the field sweep at low-temperature, the system 
is not found in the same magnetic phase as the one ob- 
tained by zero-field cooling. Such memory effect should 
be a generic feature at incommensurate to commensurate 
spin-flop transitions in spiral magnets. The magnetic ex- 
citation spectrum at the zone center in the commensu- 
rate high-field phase consists of three different modes, 
whose polarization patterns were determined by exam- 
ining different Q-positions. These results are fully con- 
sistent with the assumed field-induced flop of the spiral 
plane from be to ab plane, which explains the giant mag- 
netoelectric effect arising from the rotation of the electric 
polarization from P^ to Pa- One of the two observed c- 
polarized modes is the electromagnon of the multiferroic 
high-field phase. These results indicate that the inverse 
Dzyaloshinski-Moriya coupling, see equation (1), which 
fully explains the temperature driven multiferroic cou- 
pling, also accounts for the field-induced giant magneto- 
electric effect. 
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